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Abstract. We investigate aspects of the dynamics of a continuous atom-laser 
scheme based on the merging of independently formed atomic condensates. 
Our theoretical analysis covers the Markovian as well as the non-Markovian 
operational regimes, and is based on a semiclassical (mean-field) two-mode model. 
The role of the relative phase between the two condensates and the effect of 
interatomic interactions on the evolution of the trapped populations and the 
distribution of outcoupled atoms are discussed. 
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1. Introduction 

The dynamics of the atom laser, beyond the Born and Markov approximations have 
been investigated mostly in the framework of a rather simple model, involving a single 
condensate mode coherently coupled to a continuum of free-space modes [U O |31 IH [5J 
[6]. In recent work [7], we extended this model to a two-mode scenario. Specifically, we 
considered two initially independent Bose-Einstein condensates (BECs) consisting of 
a large number of bosonic atoms cooled into the lowest eigenmode of their respective 
traps. The two traps are assumed brought together, while atoms are outcoupled 
coherently from one of the traps only, whose condensate mode plays the role of the 
"lasing mode" which is replenished with atoms from the other trap, which serves 
as a "pump" . The chief motivation for this investigation came from the experiment 
of Chikkatur et al. [8] on the merging of BECs, which offers a conceptually simple 
scenario for a continuous atom laser. The main idea in that work relies on replenishing 
the condensate mode of the lasing trap by transfering other independently formed 
BECs using optical tweezers. 

The theoretical modelling of such a scenario has to cope with a reasonable 
description of the condensates, the non-Markovian character of the coupling to the 
atomic continuum, as well as the interatomic interactions present in the system. In an 
attempt to understand how the presence of the second condensate affects the atom- 
laser dynamics, we focused first on an interaction-free model, assuming in addition 
a vanishing relative phase between the two condensates [3. Even within these two 
simplifications, the system was found to exhibit a number of interesting features such 
as the formation of a bound state in the strong-outcoupling regime, as well as a 
dark line in the spectrum of the outcoupled atoms. In the present paper, we relax 
both of the above simplifying assumptions, addressing thus a more realistic scenario 
which takes into account inter-particlc interactions and a non- vanishing relative phase 
between the two independent BECs. As one might have anticipated, a number of new 
features emerge in the time-development of the system as well as the spectrum of the 
outcoupled atoms. 

In brief, the interplay among several parameters such as the changing chemical 
potentials due to the outcoupling, the strength of the interparticle interactions 
(nonlinearity) , as well as the strength of the outcoupling, set the stage for the 
appearance of a Josephson effect which competes with self-trapping, and also chirping, 

1. e. a variation of the energies of the outcoupled atoms. 

2. The system 

Our system consists of two independently prepared, elongated BECs (A and B), with 
N being the total number of atoms in the system. The two traps are initially far apart 
and each BEC experiences only its local potential, while only the lowest level of each 
trap (condensate mode) is populated. To allow for the merging of the two BECs, trap 
B is moved towards trap A along the radial direction x, while at the same time atoms 
are outcoupled coherently from BEC A, by applying external electromagnetic fields 
(see Ref. [7] for a detailed description of the system under consideration). Hence, to 
point out the discrete roles of the two BECs, from now on we also refer to BEC A as 
the "lasing condensate" and to BEC B as the "pumping condensate" . 

Presently, transport of BECs can be realized efhciently using optical tweezers 
which are produced by focused laser beams and offer limited trap volume and depth 
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[5]. Hence, during the merging process the two BECs can be brought as close as the 
traps' beam waist, before they start affecting each other. To be consistent with the 
experimental setup for BEC merging [S] as well as related theoretical work [§1 [10] , 
we assume two nearly identical axially symmetric harmonic traps with confining 
frequencies Uz and The atomic density profiles adiabatically follow the 

movement of the traps, while the global potential experienced by the trapped atoms 
can be modeled by a time-dependent double- well potential [7] . In general, besides the 
merging time-scale im, the details of the motion of trap B are not of great importance 
[HI [SI [in] • The crucial point is that the BEC merging must be adiabatic so that any 
kind of excitations in the system are suppressed. 

The outcoupling mechanism under consideration is based on the application of 
external electromagnetic fields which induce an atomic transition from the internal 
state of the trapped atoms to an untrapped state |/). The outcoupled atoms 
are guided through an atomic waveguide [TTl [T^] , thus resulting in an effective one- 
dimensional atom laser propagating along the weak confining axis of the waveguide. 
Such a guided atom laser has been demonstrated recently by Guerin et al. |13| and 
offers many advantages over the conventional outcoupling schemes. Formally speaking, 
the strong transverse confinement allows us to assume that the transverse dynamics of 
the free atoms adiabatically follow the slowly varying transverse potential of the optical 
guide [13] , which is assumed to be nearly the same with the transverse potential of trap 
A. In the absence of gravitational or other forces (as in the experimental setup [13]), 
the longitudinal component of the potential is practically zero and the outcoupled 
atoms propagate freely along the z direction. 



2.1. Two-mode model 

During the merging, the condensate wavefunctions start overlapping in space as the 
two traps approach each other along the x direction, and a Josephson-type tunneling 
is established between the two BECs. If the position uncertainty in the ground state of 
the traps is much smaller than the separation of the minima of the global potential, the 
overlap (and thus the Josephson coupling J) is small enough so that only the ground 
states of the traps are relevant. In first-order perturbation theory, the corresponding 
local ground-state wavefunctions iy9A(B) (r) are orthogonal and describe faithfully BEC 
A and B, at any time < t ^ [II]- If in addition the effect of interatomic 
interactions on the ground-state properties of the two wells is small i.e., if 

{ijJ^UJyLOzf/^ > Nk, (1) 

where hn is the onsite interaction energy per particle, the system under consideration 
is well-described in the framework of the standard two- mode model [H] . Introducing 
the sets of bosonic operators {a\ d} and {b\ b} for the description of the condensate 
modes A (lasing) and B (pumping) , respectively, the Hamiltonian of the system reads 

HI 

U ^ !^(a^a + b^b) + h'^UkC^Ck + hK{a)a)aa + 6^6^66) + hJ{a)b + b'^ a) 

k 

+;i^g(fc,t)(ac^ +0^4) + ^5(/c,i)(64 + Stcfe), (2) 

k k 
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where c^(cfe) is the creatioii(anmhilation) bosonic operator of a free atom with 
momentum k, and frequency 

hp 

^k = 7^, (3) 

with m the atomic mass. In deriving Hamiltonian ([2]) we have apphed the rotating- 
wave approximation for the outcouphng Hamiltonian, while we have neglected higher- 
order cross-interaction terms. Moreover, it is worth noting here that although the 
outcoupling mechanism is applied to the lasing BEC only, atoms are also outcoupled 
from the pumping condensate; albeit with a much slower rate [compare the last two 
terms of Hamiltonian ([2])]. This is due to the overlap of the two BEC wavefunctions 
and, as we will see later on, it may affect considerably the spectrum of the atom laser. 

For a weakly-interacting bosonic gas satisfying ([T]), the ground-state wavefunc- 
tions are well approximated by Gaussians and thus analytic expressions for all the 
parameters entering the Hamiltonian are readily obtained 

hott , , 

(5) 



g(fc,^) = ^V^e-^^'=/^ (6) 

where A(t) is the coupling between trapped and untrapped atomic states, r\ is the 
dimensionless separation of the traps, A — uJz/oJx, and Iz = yh/mtOz- In view of these 
relations, condition ([T]) yields the following upper bound on the total number of atoms 
we may consider in our simulations 

7V«Al/3^/2^-^, (7) 

where ott is the s-wave scattering length for collisions between trapped atoms. 

2.2. The atomic continuum 

The quadratic dependence of ujk on k is responsible for a number of mathematical 
difhculties arising in the context of atom lasers [1]. For the one-dimensional atom- 
laser model under consideration, the density of states which are available to a free 
atom diverges at frequencies close to uJc — 0, and is of the form 



where Q{uj) is the usual step function. Taking advantage of the symmetrical shape of 
the coupling and the even parity of ld^ , the spectral response of the continuum for the 
particular outcoupling mechanism under consideration reads 

sP)\(t\ 

i^H^2|g(c.,t)|VM^ ^t^^^ ^ 9H, (9) 



where g(uj,t) is readily obtained from g{k,t) using the atomic dispersion relation ([3]). 

Clearly, the spectral response ([9]) does not vary slowly for all the frequencies and 
thus we are dealing with a structured continuum which, in general, invalidates both 
Born and Markov approximations. To address fundamental mathematical difficulties 
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associated with a continuum of this type, a number of new theoretical techniques 
have been developed [H [13 HH [Uj over the last years. In the following, we adopt a 
discretization approach first developed in the context of photonic band-gap continua 
[15] , In the case of atom-lasers, the same technique has been shown capable of 
providing not only the evolution of the number of atoms in the condensates, but 
also the distribution of the outcoupled atoms in frequency domain, irrespective of the 
strength of the outcoupling and the form of the spectral response [6] . 



3. Heisenberg equations of motion 

Given the Hamiltonian ([2]) one may proceed to derive Heisenberg equations of motion 
for the operators pertaining to the two traps and the continuum: 

^ = -ii^(a) - 2iK(ataa) -iJ(6) - 2i / dkg{k,t){ck), (10a) 
dt 2 Jq 

^ = -i^(6) - 2iK{b^bb) - U{a) - 216-"' r dkg{k, t){ck), (106) 
dt 2 Jq 

%^ = -iLJk{ck)-ig{k,t){a)-ig{k,t)e-^"{b). (10c) 
dt 

Due to the presence of interatomic collisions we have the appearance of third-order 
correlation functions on the right-hand side of equations (|10aP - (|10cp . Thus, this set 
of equations is not closed, while consideration of differential equations for the third- 
order correlation functions leads to the appearance of terms of even higher order and 
so on. In general, there are no exact remedies for such mathematical problems, but an 
approximate solution can be always obtained by decorrelating higher-order correlation 
functions into products of lower ones. 

In the present work we decorrelate the third-order correlation functions appearing 
on the right-hand side of equations (|10a|) - (|10cp as follows: (a^aa) w (a^)(a)(a) and 
{b''bb) w (b''){b)(b). We thus obtain the following closed set of coupled nonhnear 
differential equations 

d(a) . fujz 

"df 

d(S) _ _ . /^^ 
dt ~ V 2 

-iLO,{ck)-ig{k,t){a)~ig{k,t)e~^-'\b). (He) 

dt 

Although such a decorrelation is, in general, a crude approximation, it turns out to 
be equivalent to a higher-order decorrelation (mean-field approximation) of the form 
(a^ad) « {a'' a) (a) and {b''bb) « when the two BECs are initially considered 

to be in coherent states. 

It is worth noting here that, in the context of the semiclassical (mean- field) theory, 
the initial coherence of the condensate(s) is preserved in time. On the contrary, 
quantum models which go beyond the mean-field approximation predict a collapse of 
the condensate wavefunction (e.g., see [131 HH])- The time of collapse tc ~ 1/(2\/]Vk) 
increases as we increase N with fixed kN, and thus the validity of the mean-field 
theory also improves. In any case, it is obvious that any experimental investigations 
pertaining to atom-laser outcoupling have to be performed for sufficiently short times 



(^ + 2«|(d)p)(a)-iJ(&)-2iy^ dkgik,t){ck), (Ha) 

l-OO 

2«|(fe)|2j (6) - i J(d) - 2ie-"' / dkg{k,t){ck). {lib) 



Effects of phase and interactions on atom-laser outcoupling from a double-well BEC6 



(i.e., for t ^ tc) so that the coherence of the atom laser is guaranteed. In this context, 
mean-field theory has been shown capable of providing rather accurate and simple 
interpretation of the experimental observations (see, for instance, pUl I^Tl [^. '23]). 

The presence of nonlinearities on the right-hand side of equations (|llap - (lllcp 
complicates somewhat their solution. Nevertheless, as long as the continuum under 
consideration is smooth (i.e., slowly varying) it can be eliminated by means of standard 
approximations (e.g., pole approximation). The atomic density of states ([8]) can be 
considered as smooth only for frequencies well above luc and for sufficiently weak 
outcoupling rates. The latter assumption, however, is not always satisfied in typical 
experimental setups (e.g., see Ref. [22l[2ll) and thus, the set of equations (|lla[) - (|llcp 
has to be treated exactly. 

To this end, we adopt here a uniform discretization approach. More precisely, 
we substitute the atomic continuum for frequencies within a range around ^2/2 (i.e., 
for < w < Wup), by a number (say M) of discrete modes, while the rest of the 
atom-mode density is treated perturbatively since it is far from resonance. The 
frequencies of the discrete modes are chosen to be tOj = js, where the mode spacing 
e is determined by the upper-limit condition of the discretization, namely Wup = Me. 
The corresponding coupling for the j mode, is determined by the spectral response ([9]) 
as follows = D{u!j)e. Discussion on the choice of Wup and the number of discrete 
modes can be found in Refs. [18] . Working as in 6^, after the discretization equations 
(fTTa| - (fTTc|) read, 

dr\ 

^ ^ -i [f,A{t)/h - S] (a) - i(J - Se-^'m - i^5.(2.), (12a) 



M 
2 ^ — A 



HBit)/h- Se-"^"^ (6) -i(J-S'e-'' )(a) -ie-" (12&) 
(cj) - i.9j (a) - i^j-e"''' (b) , (12c) 



dt 

d(g,-) 
dt 

with the shift given by S" = duj D {lj) / uj , while we have introduced the chemical 
potentials for the two BECs namely, HAit) = fujJz/2 -\- 2SK|(a(t))p and HB{t) = 
/ia'2/2-|-2fi,K|(6(t))p. In the absence of interactions, the time-dependent nonlinearities 
vanish and we recover the equations of motion for the interaction-free model of [7] 
with /iA(^) = MB(i) = fujJzl2. 

4. Simulations 

Throughout our simulations we have considered ^'^Na BECs with m — 3.818 x 10~^^ 
Kgr and att = 2.75 x 10^^ m, which are formed independently in identical harmonic 
traps with longitudinal oscillation frequency lo^ — 200 s~^ and ratio A — 0.4. 
Condition ([7]) thus yields N <C 2.5 x 10^ i.e., our two-mode model is valid for relatively 
small BECs consisting of a few hundred of atoms. 

We assume that the BECs A and B arc initially prepared in coherent states |a) 
and 1/3), respectively. Equations (|12ap - (|12el) are thus solved with initial conditions 

(a(0)) = a - V^5(0), 
(6(0))=/3=7iV/3(0)e'*(°), 

(c,(0))=0, (13) 
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where </)(0) is the initial relative phase between the two BECs. Accordingly, the initial 
number of condensed atoms in the traps A and B are given by -/Va(O) = |(a(0))p = 
Na.{Q) and iVB(O) = |(S(0))p = iV/3(0) respectively. Introducing the amplitude for 
the atomic continuum 7, at any time i > we have oi{t) + (3{t) + ^{t) = 1, so 
that iVt,.ap(t) + Nc{t) = iV, where 7Vtrap(t) = Np,{t) + N^it) and Nc{t) denote the 
population of the traps and continuum, respectively. Finally, for the sake of simplicity 
and without introducing any significant errors, the applied outcoupling pulse A(i) is 
modeled as rectangular lasting from t ~ Q to t = t. From now on, the indices A(B) 
and C, refer to traps A(B) and the continuum, respectively. 

The dynamics of the system in the framework of an interaction-free model and 
for 0(0) = were analyzed in detail elsewhere [7J. For the sake of completeness, in 
the following subsection we recapitulate briefly the main results of [7]. 

^.1. Interaction- free model and vanishing relative phase 

For the particular parameters under consideration, the weak-outcoupling regime 
corresponds to outcoupling rates A < 5 x 10^ s~^ [7] . In this regime, the system evolves 
in time with the BECs exchanging population but the oscillations are exponentially 
damped as atoms are irreversibly coupled out of the traps. The strength of the 
oscillations is determined by the Josephson coupling whereas the decay rate by the 
outcoupling strength A. Moreover, the distribution of the outcoupled atoms exhibits 
a characteristic asymmetric doublet with the two peaks separated by 2 J. 

The strong-outcoupling regime occurs for A > 5 x 10^ s~^ and is characterized 
by a continuous population exchange between the two BECs as well as between the 
BECs and the continuum. In the latter process, the pumping BEC seems to participate 
passively as its population is gradually transferred to the continuum in an irreversible 
way. On the contrary, the lasing BEC keeps exchanging population with the continuum 
even for larger times while for sufficiently strong outcoupling rates it is found only 
partially depleted in the long-time limit. The formation of such a bound state between 
the lasing condensate mode and the continuum is one of the most prominent effects 
associated with the non-Markovian nature of the dynamics. The distribution of the 
outcoupled atoms exhibits a broad peak, stemming from atoms directly outcoupled 
from the lasing BEC, while on the other hand atoms outcoupled from the pumping 
BEC give rise to a narrow peak. Most importantly, as a result of destructive quantum 
interference between the outcoupling channels in the system, the atomic distribution 
may also exhibit a dark spectral line (dip). 

We turn now to the discussion of the dynamics of the system in the case of a 
non-vanishing initial relative phase between the two BECs, as well as in the presence 
of interatomic interactions. 

4-2. Effects of relative phase 

In order to understand the role of the initial relative phase 4>{Q), we solved equations 
(|12a|) - ()12c|) numerically for increasing 0(0) € [0, 27r]. For the sake of comparison 
with the results presented in [7], throughout our simulations we have focused on an 
interaction-free model. Moreover, we have considered various outcoupling rates A, in 
order to cover both Markovian and non-Markovian dynamics. Nevertheless, before we 
proceed to such complicated simulations, it is worth recalling here that in the absence 
of atomic collisions and outcoupling (i.e., for k = and A = 0), the evolution of 
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the system is governed by the Josephson tunnchng between the two BECs only. In 
this special case, the set of equations (|12al) - (|12cp can be easily solved analytically, 
obtaining for the trapped populations 

{a^{t)a{t)) = iVA(O) cos2(Jt) + Nb{Q) sin^{Jt) + y^N^JO)N^ sin{2Jt) sin[(/)(0)] (14) 

and {b^{t)b{t)) — N — {a'^ {t)a{t)) . The first two terms in equation p^ . describe the 
population exchange between the two traps when initially iVA(O) ^ A'b(O), whereas 
for iVA(O) — A^b(O) the time evolution of the system is governed only by the third 
term which is proportional to sin[(/)(0)]. This term describes interference phenomena 
between the two BECs [3S] and was not present in the simulations performed in [7], 
as it vanishes for 0(0) = 0. In the following, we discuss the atom laser dynamics for 
various values of (t>(0). 

4-2.1. Weak outcoupling — Markovian dynamics. In the weak-outcoupling regime, 
the damped oscillatory behavior of the trapped populations reflects the Josephson 
tunneling between the two BECs, and the irreversible outcoupling of atoms into the 
continuum. As shown in [7], the evolution of the system under both of these physical 
processes strongly depends on the trap separation as well as the frequencies and the 
outcoupling rates of the two condensate modes. In addition, the present simulations 
suggest a strong dependence on the initial relative phase between the two BECs, 
which affects mostly the exchange of atoms between them. As depicted in figure [l] 
the amplitude of the oscillations varies with 0(0), while in the case of 0(0) = 7r/2, 3tt/2 
and for short times, a large fraction of the total population (more than 80%) may be 
found to occupy one of the traps only [see dashed and dotted lines in figures [11(a) and 
(b), respectively]. 

According to the inset of figure [ija), the amplitude of the oscillations varies 
sinusoidally with 0(0); a behavior which can be explained, to some extent, in terms 
of the interference term appearing in equation (|14p . which is also proportional to 
sin[0(O)]. Of course equation (|14p has been derived in the framework of a particularly 
simple model pertaining to two isolated BECs, and as such cannot describe all aspects 
of the observed dynamics. For instance, it is obvious that the amplitude variation 
depicted in the inset of figure Wi.a,) is not symmetric around 0(0) — tt while even 
the solutions we have obtained for 0(0) = and 0(0) = tt do not coincide, despite 
the fact that the interference term in equation p4|) vanishes in both cases. Such 
deviations from the perfect sinusoidal behavior can be attributed to the presence of 
the outcoupling mechanism which is not included in the simple model that equation 
(|14p refers to. More precisely, in the present scenario the two condensate modes 
are not isolated, but rather coupled to the same continuum. As a result, we have 
various outcoupling channels for the system which may also interfere, thus giving rise 
to phenomena that cannot be described in the context of isolated BECs. 

The distribution of the outcoupled atoms also varies with 0(0). More precisely, the 
characteristic doublet structure persists for any choice of 0(0), but the relative height 
of the peaks (ratio of the right to the left peak) changes periodically. For instance, 
in figure [2] we see that for 0(0) — (gray line), the right peak is more pronounced 
than the left peak, whereas the situation is the opposite for 0(0) = tt (dashed fine). 
Moreover, for 0(0) = 7r/2, 37r/2 the two peaks have nearly the same heights (solid 
black line). For a better overview of this periodic change, in the inset of figure [2] we 
present the relative height of the two peaks for various values of 0(0) (stars). The 
observed behavior can be easily understood in terms of the global symmetric and 
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antisymmetric states |±) of the double-well potential [7 . The initial populations of 



the states are P± = \y/a{0) ± W /3(0)e"*'*'' p/2 and thus their normalized ratio reads 



This expression describes the variation of the distribution of the initial population 
between the two global states with respect to 0(0). As depicted in the inset of figure 
[5J the ratio (|15p (multiplied by a constant) is in a very good agreement with the 
periodic changes observed in the relative height of the two peaks in the distribution 
of the outcoupled atoms. 

4-2.2. Strong outcoupling — Non-Markovian dynamics. Interference effects between 
the two BECs may also alter the evolution of the system in the strong outcoupling 
regime. More precisely, by varying the initial relative phase between the two BECs 
we find that the most noticeable changes in the evolution of the trapped populations 
occur for relatively short times, as well as in the long-time limit. 

Let us start with trap B for which we can identify two distinct regimes in the 
evolution of the corresponding population for 0(0) = [see gray line in figure[31[b)] . For 
short times the main part of the population is transferred into the continuum. After 
this initial transient regime, dissipation is temporarily turned off and a weak oscillatory 
population exchange between the two traps sets in. Nevertheless, trap B is empty in 
the long-time limit as its population is gradually transferred into the continuum in an 
irreversible, almost exponential, way. As depicted in figure ISjb), the initial transient 
regime gradually disappears with increasing 0(0) and for 0(0) > 57r/6, it is replaced 
by population exchange between the two BECs. In particular, for very short times the 
pumping BEC gets population from the lasing BEC, and this transfer becomes more 
pronounced as we increase 0(0) up to 47r/3. After this initial population exchange, 
the condensate mode B decays almost exponentially and is empty in the long-time 
limit. Of course, increasing further 0(0) £ (47r/3, 2tt], we find that the initial transient 
regime is gradually resettled. 

The situation is similar but more interesting in the case of trap A. More precisely, 
as depicted in figure ^a,) , the initial decay of the lasing condensate mode becomes 
faster and stronger as we increase 0(0) up to tt. For 0(0) = 0, the main part of 
the population that leaves the lasing BEC in this initial transient regime is lost into 
the continuum. The situation, however, is different for 0(0) 7^ since only part 
of the population is outcoupled while the remaining part is transferred to trap B. 
This transfer becomes more pronounced for 0(0) > 5tt/6 [see black solid and dotted 
lines in figures [3Ua,b)] . In any case, after this initial transient regime, dissipation is 
temporarily turned off and the condensate mode A exchanges population with the 
continuum only. This oscillatory population exchange persists even for longer times 
and the amplitude of the oscillations varies with 0(0). In the long-time limit, we 
have the formation of a bound state between the lasing condensate mode and the 
continuum. As a result, trap A exhibits a non-zero steady-state population which, 
as depicted in the inset of figure [3lja), also varies periodically with 0(0), and attains 
its maximum at 0(0) = tt. This behavior can be interpreted in terms of the global 
states |±) with energies ujz/i ± J. For 77 < 2.0, the Josephson coupling J > and 
thus, the antisymmetric state |— ) is closer to the edge than the symmetric state. As a 





(15) 
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result, the antisymmetric state is more protected against dissipation and contributes 
the most to the steady-state population. Moreover, as we discussed earlier, the initial 

population of the antisymmetric state is given by P_ = 0.5 — y^(5(0)/3(0) cos[0(O)] and 
attains its maximum at (/>(0) = tt. It is this periodic dependence of P- on (/)(0), that 
is also reflected in the steady-state population of trap A [see inset of figure [SJa)] ■ 

Finally, the choice of the initial relative phase may also affect significantly certain 
aspects of the distribution of outcoupled atoms in the strong-outcoupling regime. 
More precisely, the first aspect pertains to the relative height of the two peaks. As 
discussed in [7], the narrow peak stems from atoms directly outcoupled from the 
pumping BEC, whereas the broad peak is associated with atoms that emerge from 
the lasing BEC and as such experience stronger outcoupling. Thus, any variations 
in the population transfer between the two BECs with respect (/)(0), are expected to 
be reflected in the relative height of the two peaks. Indeed, as discussed above, the 
fraction of the atoms transferred from the lasing to the pumping mode increases as we 
increase 0(0) S [0, tt]. Accordingly, in this regime of relative phases the narrow peak 
related to the exponential decay of trap B becomes higher and more pronounced for 
increasing 0(0) at the detriment of the broad peak. Increasing 0(0) further, the narrow 
peak decreases again whereas the broad peak becomes stronger. The second aspect 
of the atomic distribution which is also of particular interest, involves the presence 
of a dark line (dip). As depicted in the inset of figure HJa), the dip is present only 
0(0) = or 2tt, while a new dip appears between the two peaks for 0(0) — tt. These 
observations support the explanation for the dip given in namely, it is a signature of 
destructive interference between the various outcoupling channels in the system. The 
degree as well as the nature of the interference is expected to vary with the relative 
phase between the two BECs. For instance, when A = 5 x 10^s~^, the dip appears to 
give its place to a peak which is now a signature of constructive interference between 
various outcoupling channels [see figure [3]Jb)]. 

4-3. Effects of interatomic interactions 

In general, the evolution of the system under consideration is governed by three distinct 
physical processes. More precisely, apart from the Josephson and the output coupling 
which were also present in the interaction-free model of [7] , we also have the repulsive 
collisional interactions. It is reasonable therefore to define the ratios v = kN/J and 
^ = kN/ VA which quantify the effect of interatomic interactions relative to tunneling 
and outcoupling effects, respectively. How strongly the inclusion of interactions affects 
the results obtained in the framework of the interaction-free model of 7 depends on 
these two ratios. In particular, the largest deviations are expected for > 1 and ^ > 1. 

Before proceeding with our simulations, it is worth recalling here the dynamics of 
two isolated tunnel- coupled BECs [as determined by equations (I12ap - (|12cl) for A 0]. 
In this case, the trapped population is conserved i.e., Ntia.p{t) = N, while for a 
given initial population imbalance in the system p{0) = {Na{0) — N-b{0)) / N , we have 
complete periodic oscillations between the two condensates (Josephson regime) [55] 
when [H] 

iJe(O) = '^PiOf + v/l-p(O)2cos[0(O)] < 1. (16) 

The period of the oscillations increases with increasing nonlinearity (i.e., with 
increasing N) while the oscillations become anharmonic as -ffc(O) 1. Finally, for 
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Hc{0) > 1 the oscillations are no longer complete while their period reduces with 
increasing nonlinearity {self-trapping regime) [28j . 

One can readily check that Hc{t) is a constant of motion in the absence of 
outcoupling. In our system, however, the situation is more involved as atoms are 
outcoupled from the two BECs with different rates and thus neither A'^trap nor He is 
conserved. In the following, we discuss separately the dynamics of the system for weak 
and strong outcoupling rates when it starts from the Josephson and the self-trapping 
regime. As we will see, during its evolution the system may pass from one regime to 
the other. 

4-3.1. Weak outcoupling — Markovian dynamics. In this subsection we investigate 
how the inclusion of interactions in our model affects the dynamics of the system in 
the weak-outcoupling regime. For the sake of comparison with the interaction-free 
model discussed in [7], we will focus on the case of vanishing initial relative phase 
between the two BECs. 

When the system starts from the Josephson regime [i.e., if Hc{0) < 1] and 
V < 1, the evolution of the trapped populations in the presence of interactions is 
qualitatively similar to the evolution we obtained in the context of the interaction-free 
model. The small quantitative differences depicted in figure [5] are due to interatomic 
interactions which enter the equations of motion for (a(t)) and {b{t)) as time-dependent 
nonlinearitics proportional to the corresponding trapped populations NA(t) = \ {d{t))\'^ 
and NB{t) — \(b{t))\'^ [see equations (jl2a|) and (jl2b[) ]. As a result, the chemical 
potentials of the two BECs fluctuate in time and become off-resonant. As long as < 1 
the nonlinearitics are so weak that may only shift the oscillations in the populations 
with respect to the oscillations in the case of the ideal gas, whereas for > 1 the 
nonlinearities are sufficiently strong to disturb the exchange of population between 
the two BECs. In particular, as is evident in figure [6l the oscillatory behavior of both 
occupation probabilities becomes less pronounced as we increase the nonlinearity in 
the system. Thus for sufficiently strong nonlinearities the two traps tend to decay 
nearly exponentially into the atomic continuum (see solid curves in figure [6]). 

The question that arises here is whether the system remains in the Josephson 
regime throughout its evolution or not. To answer this question, one may follow, 
for instance, the time evolution of Hc{t) [with u = KNtTap{t)/J, p{t) = {NA{t) — 
NB{t))/N, 4){t) — arg((a^6))] which is not a conserved quantity anymore, as atoms 
are continuously coupled out of the traps. In figure [7] we plot the evolution of Hc(t) 
as a function of time, for the parameters of figure [5l[a) (with N = 50). Clearly, we 
have Hc(t) < 1 for all t, and thus we may safely conclude that the system remains in 
the Josephson regime throughout its evolution. The same conclusion can be drawn 
by following the evolution of He (t) for the parameters of the other plots in figures [5] 
and [5] 

Let us consider now the situation where the system starts from the self-trapping 
regime [i.e., Hc{0) > 1] ■ As depicted in figure [51 for short-times the occupation 
probabilities for the two traps undergo partial oscillations which typically characterize 
the self-trapping regime. Gradually, however, the population imbalance between the 
two traps reduces as atoms are coupled out of trap A, e^ times faster than trap B 
and after some time we have Na = Nb (vertical gray lines). In the absence of losses 
(i.e., for A = 0) such a balanced state can be a stationary state of the nonlinear 
Schrodinger equation. In the present case, however, the system is driven again away 
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from this state as the outcoupUng mechanism is still on and the population balance is 
immediately destroyed. For sufficiently weak nonlinearities, the balanced state marks 
the transition from the self-trapping to the Josephson regime where we have complete 
damped oscillations between the two traps [e.g., see figures[5Jb,c)]. On the contrary, for 
larger nonlinearities [see figure [HJd)], the system tends to continue in the self-trapping 
regime up to the point the population imbalance between the two traps vanishes 
again. These transitions are also evident in the evolution of the corresponding He (t) , 
presented in figure [Tf^b). In contrast to figure [7ja), we see here that Hc{t) exhibits 
abrupt and instantaneous changes, which of course are not reflected in the evolution of 
the corresponding populations. However, the profile of Hc(t) (thick line) characterizes 
rather accurately the transition points discussed above. 

According to our simulations, the ratio ^ does not seem to affect considerably the 
oscillatory behavior of the occupation probabilities for the two traps. This was to be 
expected to some extent, as we are in the weak-outcoupling regime and thus these 
oscillations reflect the population exchange between the two BECs. Hence, it is the 
ratio 1/ the one that determines how strongly the population exchange is affected by 
the presence of interactions in the system. Nevertheless, our simulations show that 
the distribution of the outcoupled atoms is determined by the ratio ^. 

In figure [9] we depict the distribution of outcoupled atoms at t — 10 s for various 
atom numbers, and let us focus for a while on the upper peak of the doublet. As we 
increase the number of atoms in the system, the peak becomes broader while its shape 
starts deviating from the standard Lorentzian profile. Actually, for ^ > 1 it acquires a 
rather complicated structure. This behavior can be understood in terms of the time- 
dependent chemical potentials appearing in equations ()12ap and (I12&[) . More precisely, 
due to the applied outcoupling mechanism, the trapped population reduces with time. 
For very short times, the frequency of the outcoupled atoms that contribute to the 
upper peak is roughly /xa(0)/?i -|- J while as time goes on, yLtA(i) changes due to the 
applied outcoupling mechanism as well as the coupling to the pumping condensate. 
Hence, outcoupled atoms may involve different frequencies and the center of the atomic 
distribution is shifted accordingly towards lower or upper frequencies ( "chirping" ) ; 
a phenomenon which is reflected in the broadening of the upper peak. Moreover, 
interference effects between atoms which have the same frequency and were outcoupled 
at different times, may also give rise to unconventional atomic distributions exhibiting 
multiple spikes [see figures EJc-d)]. 

Increasing the total number of atoms in the system N or reducing the outcoupling 
rate -s/A, we essentially increase the ratio ^ and thus the effects of interactions on 
the atom-laser outcoupling are more pronounced. In this spirit one may notice, for 
example, that the same interaction effects we have just described also appear on the 
lower peak, albeit sooner. The reason is simply that this peak is associated with atoms 
outcoupled directly from trap B with outcoupling rate VXe~^ . 

In closing, we would like to point out that similar conclusions can be drawn 
by following the evolution of the atomic distribution with time, for a fixed initial 
nonlinearity i.e., for a given N. In this case, the effects of interactions become 
more prominent as time goes on, and the distribution of the outcoupled atoms 
begins exhibiting sharp spikes due to quantum interference. Moreover, all atomic 
distributions we have checked throughout our simulations show the characteristics 
discussed above, irrespective of the chosen initial conditions. Finally, it is worth 
noting that the "chirping" effect and the unconventional atom laser spectra have also 
been predicted recently by Johnsson et al. [20] , and a chirp-compensation mechanism 
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has been proposed. Moreover, chirping effects may also occur in more reahstic three- 
dimensional systems, disturbing the transverse density distribution of the atom-laser 
beam pT| . 

4-3.2. Strong outcoupling — Non-Markovian dynamics. We turn now to the discussion 
of how the inclusion of interactions in our model affects the evolution of the system 
and the distribution of outcoupled atoms in the strong-outcoupling regime. As in the 
weak-outcoupling regime, we will focus on (t>{0) = 0. 

In figures [10] and [TI] we plot the evolution of the occupation probabilities for the 
traps as functions of time. The gray curve refers to an interaction-free model whereas 
the other two curves have been obtained for a weakly interacting gas when the system 
starts from the Joscphson and the self-trapping regime. In any case, when ^ < 1 we 
may note only small deviations from the interaction- free model (see figure llOp . The 
reason is that in the strong-outcoupling regime the evolution of the system is mainly 
governed by the exchange of population between the BECs and the continuum. Thus, 
as long as ^ < 1, the effect of interactions on the evolution of the system is negligible, 
whereas when the system starts from the self-trapping regime and ^ > 1, we note large 
deviations from the interaction-free model (see black solid line in figure [TT|) . Moreover, 
for TV = 100 and the parameters of figure [TUl Hc{t) < 1 for all times, and thus the 
system remains in the Josephson regime throughout its evolution. On the contrary, for 
A'' = 200, as well as for the parameters of figure [TTl a transition from the self-trapping 
to the Josephson regime occurs rather soon i.e., at t « 9 — 13 ms. 

Our simulations also show a decrease of the steady-state population of trap 
A which is associated with the formation of a bound state between the lasing 
condensate mode and the continuum. The destruction of the bound state is due 
to the repulsive character of the interactions and has been predicted earlier by many 
authors in the context of the single- mode model [2 [5]. Here, however, we see that 
the same phenomenon also persists in the two-mode model under consideration while, 
as depicted in the inset of figure [TOT a) , the drop of the steady-state population for 
increasing atom numbers can be considered as linear to a good accuracy. 

The effect of interactions on the distribution of outcoupled atoms in the strong- 
outcoupling regime is summarized in figures 1121 For relatively weak nonlinearities 
(i.e., for ^ < 1) the distribution of the outcoupled atoms keeps all the main features 
noted in the interaction- free model [see figures [T2r a-b)] . The atomic distributions are 
mostly affected by the presence of interactions when C ^ 1 [see[T2l^c-d)]. In this case 
the chirping effect sets in and the distribution becomes broader. Moreover, quantum 
interference between atoms which have the same frequency but were outcoupled at 
different times may give rise to rich patterns involving multiple peaks and dips. 
Nevertheless, in contrast to the weak-outcoupling regime discussed earlier, these peaks 
(dips) are broader in the present case due to the strong outcoupling rates (compare 
to figure [9]). 

Finally, we would like to emphasize once more that the origin of the dips in the 
atomic distribution of figure [T^c) , is completely different from the one of the dark 
spectral line predicted in [7 . More precisely, as discussed earlier, the multiple dips 
depicted in figure 1121^ 0) stem from the interference between atoms which have the 
same energy but were outcoupled at different times (i.e., they have different phases). 
Interference patterns of this type are not expected to depend considerably on the choice 
of the initial relative phase between the two BECs, and thus they are present for any 
choice of 4>{0). On the contrary, the dark line predicted in [7] is a signature of perfectly 
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destructive quantum interference between different outcoupling channels in the system. 
This fact was also confirmed in section 14.2.21 where we showed that the dark line is 
perfect only for = 0, and 27r. Moreover, as shown in figure [T^ el. interactions also 
tend to disturb the interference phenomena between various outcoupling channels and 
thus the dark spectral line is not perfect for k 0. From the mathematical point of 
view, this is due to the time-dependent nonlinearities which render the two condensate 
modes off-resonant as well as the repulsive nature of the interactions. 

5. Conclusions 

Using a two-mode model, we have studied aspects of the dynamics of a continuous 
atom laser based on the merging of independently formed BECs. In particular, we 
examined effects of interatomic interactions and the role of the relative phase between 
the two BECs on the evolution of the trapped populations and the distribution of 
outcoupled atoms. Our simulations were performed in the weak- and the strong- 
outcoupling regimes. 

For the sake of comparison to earlier work [7] , we studied the role of the relative 
phase in the framework of an interaction-free system, and showed that it affects 
considerably the exchange of population between the two BECs, as well as the fraction 
of the population which remains trapped in the long-time limit for strong outcoupling 
rates. The distribution of outcoupled atoms exhibits always a characteristic doublet, 
with the phase determining the relative height of the two peaks. Moreover, quantum 
interference between various outcoupling channels in the system, may give rise to 
new peaks or even dark spectral lines in the distribution of outcoupled atoms. All 
of these phenomena, and in particular the unconventional form of the spectrum, 
stem from the presence of the pumping BEC in the neighborhood of the lasing 
condensate. More precisely, as long as the outcoupling mechanism is always on during 
the merging process, the transported BEC is also expected to interact with the applied 
electromagnetic radiation before the completion of the merging. 

In the absence of interactions the energies of the lasing and the pumping 
condensates are constant, whereas in the presence of interactions the two condensates 
are characterized by their chemical potentials, which depend on the number of trapped 
atoms. Hence, both of the potentials change as atoms are transferred from one 
condensate mode to the other, but mainly due to the applied outcoupling mechanism 
which affects differently the population of the two condensates. The two condensate 
modes become thus off- resonant, and this affects mostly the population exchange 
between the two modes. Of course, the most noticeable changes occur for nonlinearities 
whose strength is at least comparable to the outcoupling rate and/or the Josephson 
coupling. Moreover, for sufficiently strong interactions macroscopic quantum self- 
trapping may occur, where the population transfer from one mode to the other is 
partially turned off. This is a nonlinear effect arising from the interatomic interactions, 
and is preserved in time when the double-well system is isolated. In the present case, 
however, the presence of outcoupling affects the trapped populations and the relative 
phase between the two BECs, causing the system to enter gradually the Josephson 
regime, where it remains for the rest of its evolution. In the strong-outcoupling regime 
interatomic interactions seem to prevent the formation of a bound state between 
the lasing mode and the continuum. In particular, our simulations show that the 
steady-state trapped population decreases linearly with increasing strength of the 
nonlinearity. 
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Finally, the variation of the chemical potential of the lasing condensate with time 
is also reflected in the spectrum of the atom laser, as the outcoupled atoms experience 
a time-dependent mean-field potential, and thus their energy also varies over time 
(chirping). This effect causes a significant broadening of the distribution of outcoupled 
atoms, which may also exhibit complicated patterns reflecting the interference between 
atoms which have the same energy but were outcoupled at different time instants. 
Reducing the effect of the mean- field potential on the outcoupled atoms (e.g., by 
reducing the condensate density or working with interaction-free BEC systems in the 
vicinity of Fesbach resonances) one may in general improve the spectrum of the atom 
laser. It is worth noting, however, that an ideal continuous atoms laser working at 
steady state is not expected to suffer by chirping effects as the population of the lasing 
mode is kept constant. 
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Figure 1. Effects of the relative phase on the Markovian dynamics. The 
normahzed trapped populations are plotted as functions of time for various 
values of the initial relative phase 0(0) between the two BECs. Inset: Typical 
behavior of the amplitude of the oscillations for varying initial relative phase 
0(0). The depicted values (stars) are estimated from the evolution of the trapped 
population \{a^ (t)a{t))\'^ for various 4>{0) and in particular at its second peak. 
The corresponding amplitude for i?i(0) = is used as a reference amplitude. 
The solid and the dashed lines are sinusoidal fitting to the data, with constant 
(~ I sin[(/)(0)]|) and exponentially damped amplitude (~ | sin[?i(0)] |e-°'°5'^(°) ), 
respectively. Parameters: uiz = 200 s~^, k = 0, A = 0.4, A = 10^ s~^, and 
?7 = 1.7. Initial conditions: a(0) = 0.7, /3(0) = 0.3. Discretization parameters: 
M = 1500, ujup = 300 s-i. 




Figure 2. Effects of the relative phase on the distribution of outcoupled atoms 
in the Markovian regime. The distributions have been obtained at r = 10 s for 
various values of the initial relative phase 0(0) between the two BECs. Inset: 
Typical behavior of the relative height of the two peaks (ratio of the right to 
the left) for varying initial relative phase 0(0). The depicted values (stars) are 
estimated from the atomic distributions at t = 10 s for various 0(0). The solid 
line is a fitting to the data according to ratio (1151 . Parameters as in figure [T] 
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Figure 3. Effects of the relative phase on the non-Markovian dynamics. The 
normalized trapped populations are plotted as functions of time for various values 
of the initial relative phase (/'(O) between the two BECs. Inset: Typical behavior 
of the steady-state population of the lasing condensate mode for varying initial 
relative phase between the two BECs. The depicted values (stars) are estimated 
numerically at r = 10 s, whereas the solid line is a fitting of the form x-hy cos[<^(0)] . 
Parameters: u>z = 200 s~^, k = 0, A = 0.4, A = 4 X 10^ s"^, and r] = 1.5. Initial 
conditions: o(0) = 0.7, /3(0) = 0.3. Discretization parameters: M = 1500, 
i^up = 300 s~i. 




Figure 4. Effects of the relative phase on the distribution of outcoupled atoms 
in the non-Markovian regime. The distributions have been obtained at r = 10 s 
for various initial values of the relative phase </>(0) between the two BECs. Other 
parameters as in figure |3] but for rj = 1.7 and A = 2 X 10'^ s~^ (a); rj = 1.7 and 
A = 5 X 10^ s~'^ (b). The arrows point at the dips while the inset shows a closeup 
of the atomic distribution around it. 
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Figure 5. Effects of interatomic interactions on tiie Markovian dynamics. The 
normalized trap populations are plotted as functions of time for various atom 
numbers and uJz = 200 s~^, A = 0.4, A = 10^ s~^, and ri = 1.5. The gray curves 
correspond to an interaction-free Bose gas. Initial conditions: Josephson regime 
with a{0) = 0.7, /3(0) = 0.3, </>(0) = 0. Discretization parameters: M = 1500, 
i^up = 300 . 
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Figure 6. As in figure [5l for rj = 1.7. 
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Figure 7. The evolution of He (thin Une) and its envelope (thick line), are plotted 
as functions of time for the parameters of figure [5ja) {N = 50) and figure [8ld). 
The arrows point at the regions where the system passes from one regime to the 
other. 
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Figure 8. Effects of interatomic interactions on the Markovian dynamics. The 
normalized trap populations are plotted as functions of time for various atom 
numbers: (a) N = 50; (b) N = 100; (c) TV = 150; (d) N = 200. The vertical gray 
line denotes the border between self-trapping and Josephson regime. Parameters: 
ijJz = 200 s~^, A = 0.4, A = 10^ s~^, = 1.7. Initial conditions: self-trapping 
regime with o(0) = 0.7, /3(0) = 0.3, </'(0) = 0. Discretization parameters: 
M = 1500, cjup = 300 s-i. 
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Figure 9. Eflfects of interatomic interactions on the distribution of outcoupled 
atoms in the Markovian regime. The distribution (a) corresponds to an 
interaction-free Bose gas whereas the other three distributions have been obtained 
for a weaJjly interacting Bose gas with AT = 25 (b), A' = 50 (c), AT = 100 
(d). Parameters: uiz = 200 s-\ A = 0.4, A = 10^ s-f, ?? = 1.5, r = 10 s. 
Initial conditions: Josephson regime with a(0) = 0.7, ,9(0) = 0.3, (j>{0) = 0. 
Discretization parameters: M = 1500, cdup = 300 s~^. 
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Figure 10. Effects of interatomic interactions on the non-Markovian dynamics. 
The normalized trap populations are plotted as functions of time for u)^ = 200 s~^, 
A = 0.4, A = 4 X 10'^ s~^, ri = 1.5. The inset shows the typical behavior of the 
steady-state population of trap A at r = 10 s for various atom numbers. The gray 
curves and circles correspond to an interaction-free Bose gas. Initial conditions: 
Josephson regime {N = 100) and self-trapping regime (A'^ = 200) with a(0) = 0.7, 
/3(0) = 0.3, 0(0) = 0. Discretization parameters: M = 1500, oJup = 300 s~^. 
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Figure 11. As in figure [Tol starting from the self-trapping regime for A = 
2 X 10^ s"2 and T) = 1.6. 




Figure 12. Effects of interatomic interactions on the distribution of outcoupled 
atoms in the non-Markovian regime, (a-d) Atomic distributions at r = 10 s for 
weakly interacting Bose gas with Af = 50 (a), N = 100 (b), N = 200 (c) and 
N = 300 (d). Parameters: lj^ = 200 s"!, A = 0.4, A = 2 X 10^ s-2, = 1.7. (e) 
A closeup of the atomic distribution around the dip for A = 10'^ s~^. The gray 
curves correspond to an interaction-free Bose gas. Initial conditions: a(0) = 0.7, 
/3(0) = 0.3, 0(0) = 0. Discretization parameters: M = 1500, Wnp = 300 s'^. 



